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In this letter we study the fuU semi-conservative treatment of a model for the co-evolution of a virus 
and an adaptive immune system. Regions of viability are calculated for both conservatively and semi- 
conservatively replicating viruses interacting with a realistic semi-conservatively replicating immune 
system. The conservative virus is found to have a selective advantage in the form of an ability to 
survive in regions with a wider range of mutation rates than its semi-conservative counterpart. This 
may help explain the existence of a rich range of viruses with conservatively replicating genomes, a 
trait which is found nowhere else in nature. 
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DNA is often called the molecule of life. The vast ma- 
jority of organisms in nature store their genetic informa- 
tion in the form of double stranded DNA, which provides 
a number of benefits over its close relative and likely pre- 
decessor, RNA^ . These benefits include a predictable sec- 
ondary structure and a resistance to auto-catalytic cleav- 
age and hence a longer half-life. However, viruses stand 
out as a notable exception to the "DNA world" , employ- 
ing a variety of methods for genetic storage including 
single stranded RNA (e.g. tobacco mosaic virus), dou- 
ble stranded RNA (such as reovirus), linear and circular 
single stranded DNA (including parvovirus and bacte- 
riophage (/)X174, respectively) and a variety of double 
stranded DNA types (examples of which include bacte- 
riophage T4, polyoma virus and poxvirus). While dou- 
ble stranded DNA replicates semi-conservatively, many 
of these viruses replicate conservatively, wherein multi- 
ple copies of a single strand are made, and the original 
strand is conserved. This suggests that either the bene- 
fits of DNA must be less for viral species (possibly due 
to a shorter genome length or unusual life cycles) and/or 
a selective advantage must exist for conservatively repli- 
cating viruses. In this paper we use the the quasispecies 
model to search for such an advantage. 

Introduced by Eigen in 1971, the quasispecies modelSi^ 
has been used to study various characteristics of conser- 
vatively replicating organisms ranging from equilibrium 
data to punctuated evolution'^i^i^i'^i^i^i^°i"i^^i^3i^'^i^^ The 
model consists of a population of independently replicat- 
ing genomes (/>, each of which is made up of a set of 
"letters" siS2---Sn chosen from an "alphabet" of size 
S. S is usually chosen to be two for simplicity or four, 
as in this paper, to model the nucleotides A, C, T and 
G. Each possible genome is assigned a fitness that dic- 
tates its fecundity. This mapping of fitness to geno- 
type can be represented by a unique "fitness landscape" . 
The process of replication includes a probability of point 
mutation per base pair e that is generally assumed to 
be genome-independent. By associating phenotype with 
genotype and assuming a first-order dependence of the 
growth rate on concentration, a set of differential equa- 
tions can be solved to describe the competition between 
various genotypes^. Although the model incorporates 



numerous approximations, it is well suited to describing 
small RNA genomes and viruses and many of its pre- 
dictions have been experimentally verified. One of the 
major successes of the model lies in recent work on novel 
anti- viral therapiesiS^. 

The quasispecies model has recently been extended to 
the coevolution of hosts and parasites and the particu- 
lar case of an adaptive immune system interacting with a 
virusiS^. Viruses make detrimental use of host biochem- 
ical processes while the immune system expends enor- 
mous effort to keep viral concentrations as low as possi- 
ble. As the immune system develops new defenses, the 
virus must adapt to defeat them. The immune system 
must then evolve to destroy the newly resistant strains, 
and a non-linear co-evolving feedback loop is created. 
To model this behavior, the immune system and virus 
are both assumed to evolve on a single fitness peak land- 
scape, where the fitness of all genomes are equal with the 
sole exception of a single master sequence of far greater 
fitness, or 



A(0) = 



(1) 



where A{(l)) represents the fitness of genome 4> ^^^d 
{cr, 77} = {dv^riv} for the virus and {cr, r?} ee {ais,r]is} 
for the immune system. This landscape is dynamic in 
that the fitness peak is allowed to shift from one genome 
to another at specified intervals. To model the inter- 
species interaction, the master sequence for the immune 
system is assumed to impose a death rate 5 on the cor- 
responding viral sequence. When this coincides with the 
viral master sequence, the viral fitness peak shifts out of 
self-preservation. The new viral master sequence regen- 
erates on a time scale , defined as the time required for 
the new master sequence to outnumber the old. At this 
point, the immune fitness peak shifts to match the vi- 
ral master sequence and will regenerate on a similarly 
defined timescale Tis after which the viral peak shifts 
again. These steps are iterated so that the virus traverses 
genome space with the immune system in hot pursuit. 
Using recent results on the dynamics of a quasispecies on 
time-dependent landscapes^, Kamp and BornholdliiSiiS 
found expressions for the long-term survival of a conser- 
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vatively replicating virus and immune systemi^iiS by con- 
sidering the behavior of each species on a dynamic fitness 
landscape. In essence, every time the fitness peak shifts, 
the concentration of the master sequence drops dramati- 
cally and begins to regrow, while the large concentration 
of the old master sequence drains away. Rigorously, the 
dynamics of a quasispecies on a single fitness peak land- 
scape follow a set of differential equations defining the 
evolution of the various genomes in terms of their Ham- 
ming distance HD{(j), (/>'). This is defined as the smallest 
number of point mutations that lead from to 0'. Uti- 
lizing the radial symmetry of the fitness landscape, the 
quasispecies equations become 
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where x G {u, is}, wi is the concentration of sequences of 
Hamming distance / from the master sequence, fis{t) = 
Y,i Aisil)wi^is = {(7is-r]is)wo^is+'r]is for the immune sys- 
tem, fvit) = S for the viral sequence that coincides with 
the immune master sequence and fvit) — otherwise, e 
represents the point mutation probability (more complex 
mutations such as insertions and deletions, as well as the 
possibility of recombination, are ignored). 

Building on recent work regarding the dynamics of 
a quasispecies on mobile fitness landscapes^, a num- 
ber of qualitative features have come to light. Besides 
the standard (albeit modified by the dynamic nature of 
the landscape) error catastrophe at high e^,, the virus is 
only viable above a given minimum mutation rate, below 
which it is unable to keep up with the moving landscape 
(dubbed the "adaptability" catastrophe). At these low 
mutation rates, each shift of the landscape is followed 
by a period of time wherein the new master sequence 
attempts to build up an equilibrium distribution. How- 
ever, before the new master sequence can rebuild to the 
levels reached by the old master before the shift, the 
fitness peak moves again. As this process repeats, any 
initial quasispecies will disappear and genomes become 
stochastically distributed. Further, an optimal mutation 
rate for the immune system can be found that minimizes 
the range of viral mutation rates that allow persistence of 
the viral quasispecies. This optimal rate has been deter- 
mined and found to be independent of the parameters of 
the model and the properties of the viral species as well 
as comparing admirably with the rates of somatic hyper- 
mutation in B-cells^^. The agreement suggests that this 
model, although approximate in nature, captures the ro- 
bust features of the co-evolution of viruses with an adap- 
tive immune system. 



FIG. 1: A schematic model of DNA replication. Adapted 
fromSa. 



Although the host-parasite model has yielded impres- 
sive successes, work has been restricted to conserva- 
tively replicating systems, which differ greatly from the 
true semi-conservative systems that dominate nature. In 
a conservative system, multiple, possibly error prone, 
copies of an original strand are produced without harm- 
ing or changing the original strand. Thus, the original 
quasispecies model is ideally suited for the study of RNA 
viruses or in-vitro RNA evolution experiments. Semi- 
conservative replication follows a different route shown 
schematically in Fig. 1. DNA exists as a tightly bound 
double helix structure, where each strand is connected 
to a complementary strand <p, where (/> represents the 
complement of strand ^, and can be written as Sj^ S2 ■ ■ ■ s„ 
where we assign the nucleotides A = 1,G = 2,T=3,C = 
4 and Sj — (st + S/2) mod S. In order to undergo repli- 
cation, the double helix unzips to free two single strands 
(j) and (j). Each of these is replicated to produce an error- 
prone complement, yielding {0, </)'} and {0,0'}. Each 
error must result in a base mismatch, which can be rec- 
ognized and selectively repaired by enzymes in the cell. 
These enzymes can recognize the new strand and ensure 
that the mismatch is repaired by replacing the new mis- 
matched base, keeping the effective error rate e low. In 
the final stage, the strands become indistinguishable and 
various maintenance enzymes repair the remaining errors 
with a 50% probability of correcting the mismatch in ei- 
ther strand. Thus, the final result are two new pairs, each 
consisting of two new strands, {0",0"} and {(/)"',(/)"'}. 
Recent work by Tannenbaum et. aS^ has extended the 
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quasispecies model to the case of semi-conservative repli- 
cation, which was found to display significantly different 
behavior in the infinite time limit on a static landscape. 

To properly treat a semi-conservative quasispecies 
model on a single fitness landscape, ignoring back mu- 
tations, Eqn. (2) must be recast as^ 
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where/ (t) is defined above. Here, we examine the dynam- 
ics of the semi-conservative equations within the confines 
of Kamp and Bornholdt's model of co-evolution. This 
study will focus on the dynamical aspects of Eqn. (3), 
as opposed to the equilibrium effects studied by Tannen- 
baum et. al.^^ . Following Kamp and Bornholdt, the 
condition for the viability of the quasispecies is 



> h 



(4) 



where k represents the ratio of master sequence concen- 
trations at the beginning and end of an entire cycle of 
landscape shifts in an unconstrained system compared 
to the equivalent growth of a random sequence far from 
the peak. Obviously, if the master sequence outgrows the 
random sequence over this period, even with the concen- 
tration losses incurred by the peak shift, it will survive 
for all times. If the master sequence is outgrown by the 
random sequence, i.e. if k < 1, the master sequence will 
diminish and disappear at long times. 

After a fair bit of work, the condition for viability of 
the immune genome can be expressed as^^ 

fT,sQ.(l-Q./2)""-i 
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A conservatively replicating virus interacting with a 
semi-conservative immune system will follow the behav- 
ior described by the conservative model of Kamp and 
Bornholdt, 



{S - l){ay - r]y)qy 



, (10) 



where q = 1 — e represents the replicative fidelty per base 
pair, T — Tis + Ty as always, and 



(11) 



In contrast to previous work, is described by the semi- 
conservative result described by Eqn. (7). In Fig. 2 we 
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FIG. 2: Regions of viability for a co-evolving host-parasite 
system. The contours for both semi-conservative and conser- 
vative virus interacting with a conservative immune system 
are shown. The vertical lines represent the region where only 
the conservative virus is viable, while the horizontal lines rep- 
resent the region where the immune system is only viable for 
a conservative virus, riis — Uv — 50, ais = av = 100, rjis = 
nv = 1,5 = 200. 

plot the regions of viability defined by Ky^is > 1 with a 
particular set of parameters for a semi-conservative im- 
mune receptor interacting with both a conservative and 
semi-conservative virus. It is immediately clear that the 
conservative virus can survive under a wider range of 
conditions than the semi-conservative, thus obtaining a 
selective evolutionary advantage. While the immune sys- 
tem is slightly less robust for the semi-conservative virus, 
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this effect is small and lies near the region where the im- 
mune system becomes unviable, and hence should have 
little effect on real systems. Although this behavior is 
dependent on the parameters of the model, the qualita- 
tive trend was robust for the vast majority of biologically 
reasonable parameter choices. Further, the increased vi- 
ability lies both near and far from the optimal immune 
mutation rate, important since a population of viruses in 
nature is expected to interact with a range of immune 
system properties, and the behavior away from the op- 
timum should play an important role in the evolution of 
the system. Thus, we have demonstrated the existence 
of a clear selective advantage for conservative replication 
in viral species. 

We conclude with a cautionary note. One must al- 



ways take great care in extracting grand predictions from 
simplified models such as this one. In nature, a myriad 

of evolutionary pressures battle for dominance and it is 
often difficult to pinpoint the selective advantage for a 
given trait. However, the model presented here likely 
captures many robust features of the evolution, and sug- 
gests a possible explanation for the success of conserva- 
tive viruses in nature. 



Acknowledgments 

The authors gratefully acknowledge useful discussions 
with Emmanuel Tannenbaum and Brian Dominy. 



^ D. Voet and J. Voet, Biochemistry (John Wiley and Sons, 

Inc., New York, NY, 1995). 
^ M. Eigen, Naturwissenschaften 58, 465 (1971). 
^ M. Eigen, J. McCaskill, and P. Schuster, Adv. Chem. Phys. 

75, 149 (1989). 

* M. Nilsson and N. Snoad, Phys. Rev. Lett. 84, 191 (2000). 

B. Drosscl, Adv. in Phys. 50, 209 (2001). 

C. O. Wilke, C. Ronnewinkel, and T. Martinetz, Phys. 
Rep. 349, 395 (2001). 

^ P. Tarazona, Phys. Rev. A 45, 6038 (1992). 

* S. Galluccio, Phys. Rev. E 56, 4526 (1997). 

^ P. Schuster and P. F. Stadler, Computers Chem. 18, 295 
(1994). 

° L. Peliti, cond-mat/9712027 (1997). 

^ S. Altmeyer and J. McCaskill, Phys. Rev. Lett. 86, 5819 
(2001). 

^ P. R. A. Campos and J. F. Fontanari, Phys. Rev. E 58, 

2664 (1998). 



D. Alvcs and ,]. F. Fontanari, Phys. Rev. E 57, 7008 (1998). 
C. O. Wilkc, C. Ronnewinkel, and T. Martinetz, Phys. 
Rep. 349, 395 (2001). 

J. Krug and C. Karl, Physica A 318, 137 (2003). 

S. Crotty, C. E. Cameron, and R. Andino, Proc. Nat. Acad. 

Sci. 98, 6895 (2001). 

L. A. Loeb, J. M. Essigmann, F. Kazazi, J. Zhang, K. D. 
Rose, and J. I. MulUns, Proc. Nat. Acad. Sci. 96, 1492 
(1999). 

C. Kamp and S. Bornholdt, Phys. Rev. Lett. 88, 068104 

(2002). 

C. Kamp, C. O. Wilke, C. Adami, and S. Bornholdt, Com- 
plexity 8, 28 (2003). 

E. Tannenbaum, E. J. Deeds, and E. I. Shakhnovich, Phys. 
Rev. E (2004), in press. 

Y. Brumer and E. Shakhnovich (2004), to be published. 



